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ABSTRACT 

We present an analysis of the asymmetries in the population of Galactic M-giant stars 
present in the 2MASS All Sky catalogue. Several large-scale asymmetries are detected, 
the most significant of which is a strong elliptical-shaped stellar over-density, close to 
the Galactic plane at {£ = 240°, b = —8°), in the constellation of Canis Major. A small 
grouping of globular clusters (NGC 1851, NGC 1904, NGC 2298, and NGC 2808), 
coincident in position and radial velocity, surround this structure, as do a number of 
open clusters. The population of M-giant stars in this over-density is similar in number 
to that in the core of the Sagittarius dwarf galaxy. We argue that this object is the 
likely dwarf galaxy progenitor of the ring-like structure that has recently been found at 
the edge of the Galactic disk. A numerical study of the tidal disruption of an accreted 
dwarf galaxy is presented. The simulated debris fits well the extant position, distance 
and velocity information on the "Galactic Ring" , as well as that of the M-giant over- 
densities, suggesting that all these structures are the consequence of a single accretion 
event. The disrupted dwarf galaxy stream orbits close to the Galactic Plane, with a 
pericentre at approximately the Solar circle, an orbital eccentricity similar to that of 
stars in the Galactic thick disk, as well as a vertical scale height similar to that of 
the thick disk. This finding strongly suggests that the Canis Major dwarf galaxy is a 
building block of the Galactic thick disk, that the thick disk is continually growing, 
even up to the present time, and that thick disk globular clusters were accreted onto 
the Milky Way from dwarf galaxies in co-planar orbits. 
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1 INTRODUCTION 

It is now generally accepted that galaxies assembled by 
the h i erarchical merging o f dark matter halos f White, Reea 
Il978l: IWhite. Frenklll99lh . with the accretion of luminous 
(dwarf) galaxies playing a significant, or probably even a 
major, role in the formation of the visible structures. In 
the Milky Way, this picture is dramatically illustrated by 
the existence of the Sagittarius dwarf spheroidal galaxy, 
which is in an advanced state of tidal disrupti on, and has 
popu lated the outer halo with M gian t stars dlbata et alJ 
l2002l). carbon stars ilbat a et al.|r 2001al'l and globular clus- 
ters bellazzini. Ferraro i^^8^ir2003l) . 

The accretion of the Sagittarius galaxy is currently 
the only strong evidence that the Milky Way is absorb- 
ing satellite galaxies. However, the ring-like structure that 
has recently been discovered around the Galaxy could be 
the consequence of another such event. First discovered 



as an o ver-density in blue st ars of the Sloan Digital Sky 
Survey l)Newberg et alJl2nn2h . this "Ring" has since been 
probed using phot ometry from the INT Wide Field Camera 
ilbata et al.l 20 03|. hereafter 10 3). SDSS spectrometry and 
Dhotomet rv (lYannv et al.ll2n0.1 he reafter Y03) and 2MASS 
M-giants l|Fi,ocha-Pinto et al.ll2003( ). Surrounding the Galac- 
tic disk, its Galactocentric distance ranges from ~ 15 kpc to 
~ 20 kpc in fields taken within 30° of the Galactic plane and 
for 122° <£< 227°. 

It has been proposed (103 and Y03) that the structure 
could be the tidal stream stripped away by the Milky Way 
tides from a satellite galaxy that has an orbital plane close 
to the plane of the Galaxy. In light of thi s interpretation, it is 
interesting to reexamine the lAbadi et a l. (2002) simulation 
of the formation of a disk galaxy assembled hierarchically in 
a ACDM cosmology. They showed that the thin and thick 
disks of their simulated galaxy, which were similar to those 
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of the Milky Way, do not share the same origin. The simu- 
lated thick disk component was constructed mainly through 
accretions of satellites with an orbit close to the Galactic 
plane. These satellites are dense enough for their orbits to 
circ ularize before they are disrupted. 

iHelmi et alJ (|200,?ll also used this simulation to follow 
more closely the behavior of an in-plane accretion. They 
concluded that it naturally leads to the formation of a ring- 
lik e structure around th e galaxy similar to the one found 
by iNewberg et sl\ l)2002h . with two possible explanations. 
One possibility is that the ring is a tidal arc produced by 
stars stripped away from the parent satellite during a re- 
cent pericentric passage. This would produce an asymmetric 
structure above and below the disk, limited in Galactic lon- 
gitude and with a significant velocity gradient in the Galac- 
tic longitude direction. Alternatively, it could be a shell-like 
structure, produced by ancient minor mergers, similar to the 
shells observed in elliptical galaxies. In this case, the ring 
would be rather symmetric and with no velocity gradient. 

To discriminate between these competing scenarios, it is 
essential to analyze a large-scale, homogeneous data set. Un- 
til now, this was the main limitation, but with the release of 
the 2 Micron All Sky Survey, All Sky Release (2MASS ASR) 
catalogue we now have an invaluable tool to determine the 
characteristics of this Galactic "Ring" . Apart from covering 
all the sky, the 2MASS ASR contains precise near infrared 
photometry, so it is not much hampered by extinction in the 
low Galactic latitude regions with high dust-contamination 
where the "Ring" is primarily located. 

We use this catalogue to probe the ring-like structure 
and determine its dimensions and characteristics on the ce- 
lestial hemisphere opposed to the Galactic centre. In sec- 
tion 2, we discuss the data set we used for our study, pre- 
sented in section 3. In section 4, we present the results of 
our search of clusters in the ring-like structure. Section 5 
presents N-body simulations of the tidal disruption of the 
progenitor of the observed structures. Finally, in section 6, 
we discuss our results and the consequences on the origin 
of the "Ring". Throughout this work, we assume that the 
Solar radius is Rq = 8kpc, that the LSR circular veloc- 
ity is 220kms~^, and that the peculiar motion of the Sun 
is (Up = 10.00 kms-\Vo = 5.25kms-\W^o = 7.17kms-^ 
iDehnen fc Binnevlll998bll . 



2 DATA 

The low Galactic latitude regions that we study here are 
mainly hampered by two effects: by reddening due to the 
high concentration of dust in the Galactic plane and by the 
presence of numerous foreground disk stars. The 2MASS 
near infrared photometry has the ability to probe the highly 
obscured regions, but it cannot be directly used to remove 
the foreground disk stars. Therefore, in order to diminish the 
contamination of the closest disk stars, we cross-identified 
the 2MASS and the USNO-Bl.O catalogues in order to com- 
plement the 2MASS photometry with proper motion data. 
This was done for stars with |6| > 5°. 

In the present study we will work with two 
samples of stars: sample A is simply the 2MASS 
ASR point-source catalogue (with 2MASS flags 
bl_flg = 111, cc_flg = 000, gal.contam = 0,mp_flg = 0); 



It will be used to search for the Galactic "Ring" population. 
Sample B is identical to sample A, except that we have 
removed all sources with significant {3a) proper motions. 
The contamination fraction of objects with proper motion 
corresponds to ^ 35% of the 2MASS sources, but since they 
are irregularly distributed in a volume of a few kiloparsecs 
around the Sun this sample cannot be used to detect 
assymetries. The proper motion discrimination is however 
very useful for plotting colour-magnitude diagrams: the disk 
dwarfs are particularly affected by this process and, since 
they tend to blend with red giant populations in CMDs, 
filtering them allows us to follow distant red giant branch 
(RGB) populations up to 2 magnitudes deeper. However, 
RGB stars beyond ~ 5 kpc are not removed with th is filter, 
judgin g from our tests in globular cluster fields iMartinI 

We de-redden the photome tric data using the 
ISchleeel. Finkbeiner fc Davij il998l) algorithm which inter- 
polates the E(B — V) value from the IRAS 100 micron emis- 
sion all-sky maps. All the J,H,Ks magnitudes stated below 
are extinction-corrected. In order to use only precise pho- 
tometry, we remove sources with Kg > 14.3 (the 10-sigma 
limit of the 2MASS photometry) and with E(B - V) > 0.555 
(which corresponds to E(J — Kg) > 0.3). This last limit 
mainly removes zones near the Galactic plane with |&| < 5° 
from our samples. To make sure we can compare quanti- 
tatively the distribution of sources in both hemispheres, we 
divide the sky into 0.1° x0.1° pixels (in Galactic coordinates) 
fiagging as unusable each zone in which the extinction ex- 
ceeds the chosen limit. However, to maintain North-South 
symmetry we also fiag as unusable the symmetric pixel on 
the opposite side of the Galactic plane. We also limit our 
present samples to \b\ < 40°. 



3 RESULTS 

3.1 M-giant sample 

From sample A, we first select a subsa mple of late M- 
giant stars using the same colour limits as lMaiewski et all 
f2003): 0.85 < J - < 1.3, J - H < 0.561(J - K^) + 0.36 
and J — H > 0.561(J — Ks) -I- 0.22. We measure the distance 
of the Galactic M-giants using the RGB o f the Sagittarius 
dwarf galaxy as a reference, as calibrated bv'Maiews ki et all 
12003): Ks = -8.650(J - Ks)o + 20.374. By assuming a dis- 
tance modulus of m — M = 16.9 for th e stellar populatio n 
at the core of the Sagittarius galaxy jlbata et al. lll99^ . 
one can estimate the distance to other RGB populations. 
Figure 1 compares the 2MASS colour-magnitude diagrams 
(CMDs) of the Sagittarius dwarf with those of the Large 
Magellanic Cloud (LMC) and the globular cluster Omega 
Centauri; the corresponding distance calibrations are also 
displayed. The relatively poor fits for the LMC and Omega 
Cen, are probably due to differences in age and metallic- 
ity between these objects and the Sgr galaxy. Using the Sgr 
RGB calibration to measure the distances of the LMC and 
Omega Cen would incur a ~ 30% systematic error, with 
these distances being underestimated. This systematic un- 
certainty will cloud all of the distances derived below for the 
Galactic M-giant population. 

A comparison between the northern and southern hemi- 
spheres of the distribution of source distances, estimated 
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Figure 1. Comparison of colour-magnitude diagrams of our M- 
giant distance calibrator, the Sagittarius dwarf galaxy (left-hand 
panel), with the LMC (middle panel) and Omega Centauri (right- 
hand panel). The proper-motion cleaned catalogue (sample B) is 
used here. The Sgr CMD is derived from a 1° circle, centred on the 
core of the galaxy, the LMC CMD comes from an annulus between 
5° and 6° , while the Omega Cen CMD is from data within 20' of 
the core of the cluster. The J — Kg, Kg colour- magnitude selection 
window is shown between the dashed vertical lines (though note 
the additional selection in J — H colour is not displayed here). The 
arrow shows the reddening vector, with a length corresponding to 
the chosen reddening limit of E(B — V) < 0.555. In th e left-hand 
panel, the solid line shows the lMaiewski et al.| RGB cali- 

bration for the Sgr dwarf (m - M = 16.9: llbata et al. Il997l) . The 
parallel sequence approximately two magnitudes brighter than 
the Sgr RGB is the RGB of the bulge. The solid line in the middle 
and right-hand panels show the corresponding ridge line assum- 
ing a distance modulus of 18.51 5 for the LM C IClementini et aD 
l2n0.?r and 13.53 for Omega Cen <Harrisll99fl) . Using the Sgr fidu- 
cial to determine the distance to the LMC or Omega Cen would 
lead to a ~ 30% underestimate of their distances. The substantial 
population red- ward of J — Kg = 1.3 in the LMC are AGB carbon 
stars, which highlights the need for the upper J — Kg colour limit. 



from sample A, using the Sgr RGB fiducial described above, 
shows two striking asymmetries, one in each hemisphere 
(see Figure 2). The southern over-density, located between 
220° < I < 260° is the most significant (detected with 
S/N = 30). It ranges from 10.0 kpc to 15.8 kpc from the 
Galactic centre and is centred at ^ 12.6 kpc. The north- 
ern over-density ranging from I = 140° to I = 220° is fur- 
ther away, between 15.8 kpc and 20.9 kpc, with its maxi- 
mum between 16.6 kpc and 19.0 kpc (this peak is detected 
at S/N = 13). In both cases there are approximatively twice 
as many stars in the over-density than in the other hemi- 
sphere. 

To examine the structures creating these two over- 
densities, the data in sample A was binned in Galactic co- 
ordinates, as presented in Figure 3. In those pixels where 
more than 50% of the area is below the extinction limit, 
we have corrected the counts for the area incompleteness; 
otherwise, the pixel values are set to zero. We impose 
North-South symmetry to be able to compare the counts 
between the hemispheres. In the southern hemisphere, for 
the 15.2 < (m — M)gc < 15.7 range, an object is centred at 
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Figure 2. The Galactocentric distance modulus (m — M)gc 
distribution for M-giant stars in sample A with |6| < 40°, 
E(B — V) < 0.555 and for different ranges in t. Note that in 
this, and all subsequent figures, those regions removed by the 
extinction cut have their corresponding symmetric region about 
the Galactic plane removed as well. The star distribution in 
the southern hemisphere (full line) is compared to the north- 
ern hemisphere distribution (dashed line); the difference between 
these distributions is also plotted (dotted line). The 220° < 
I < 260° region (left panel) shows an over-density in the south- 
ern hemisphere, ranging over 15.0 < (m — M)qc < 16.0 with a 
maximum at (m — M)qc = t5.5 (the signal to noise ratio of 
this peak is S/N = 30). The 140° < £ < 220° region (right 
panel) shows an over-density in the northern hemisphere, rang- 
ing over 16.0 < (m — M)qc < 16.6 with a maximum between 
(m - M)gc = 16.1 and (m - M)gc = 16.4 (this peak has S/N = 
13). 



£ 240° and b ~ —8°, in the direction of the Canis Major 
constellation. It is very dense at its core — with up to eight 
times the number of stars compared to the symmetric zone 
in the northern hemisphere — and it seems to be approxi- 
mately elliptical. These structures are more clearly seen in 
the difference image displayed in Figure 4, where we have 
subtracted the star-counts from the opposite hemisphere. 
The Canis Major over-density is striking, as is the wide arc- 
like structure in the northern hemisphere. The location of 
these features is made clearer with schematic overlays in 
Figure 5. 

To measure the dimensions of the Canis Major ob- 
ject, we examine the profiles in the latitude and Heliocen- 
tric distance directions (Figure 6). Analysis of the verti- 
cal distribution of the population is hampered by the lack 
of non-dust-contaminated regions for |fe| < 7°. This puts 
a high degree of uncertainty on the profile fits since the 
population could correspond to a Gaussian distribution cen- 
tred at 6 = -6.3° ± 1.0° with a FWHM vertical height of 
1.6 ± 0.2 kpc, or alternatively, an exponential distribution 
with a 0.73 ± 0.05 kpc scale height. It is interesting to note 
that this scale height derived for the Canis Major structure 
at £ — 240° is statistically identical to that obtained by 103 
for the Galactic Ring, in a field at ^ = 123°. This lends 
some support to the possibility that the structures are re- 
lated. The Heliocentric distribution of the object can be fit 
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Figure 3. The distribution of M-giant stars (from sample A) for 
different ranges in Galactocentric distance modulus (m — M)gc- 
(We choose to display these distributions in Galactocentric coor- 
dinates because the over-densities are more peaked than in He- 
liocentric coordinates, whereas the use of the distance modulus 
makes it easier to visualize the effects of extinction and popu- 
lation differences). In the upper two panels a pixel corresponds 
to a 4° X 2° rectangle, while in the bottom panel a pixel covers 
8° X 4°; the scale colour is different for each panel to emphasize 
the structures. The circular hole at ^ ~ 280° and b ~ —35° comes 
from the removal of the Large Magellanic Cloud; due to the sym- 
metry imposed on the map, this hole is reflected in the North 
as well. The top panel contains the sources of our sample with 
15.2 < (m — M)gc < 15.7. The most significant over-density in 
this distance range is created by a dense object at 220° ^ ^ ^ 
260° and h > -20°, towards Canis Major. The middle panel 
shows the distribution of sources with 15.7 < (m — M)qc < 16.1, 
while the bottom panel shows 16.1 < (m — M)gc < 16.7. In the 
top, middle and bottom panels, a black pixel corresponds to 40, 
16 and 16 counts, respectively. 

by a Gaussian centred at = 7.1 ± 0.1 kpc with FWHM 
of 4.2 ± 0.3 kpc. While the centre of the distribution can 
be measured reliably, its FWHM should be taken with cau- 
tion since the distribution is skewed to large distances. This 
could be due to the angle of observation of the structure if 
it is elongated in a direction that is not perpendicular to the 
line of sight. The precise extent of this object along the lon- 
gitude axis is harder to compute because the observations 
arc hampered by high extinction regions and by the rapid 
increase in the number of disk stars at low latitude. 

The Northern Arc over-density is visible for 
16.1 < (m- M)gc < 16.7 (bottom panel of Figures 3 
and 4). It begins near the Galactic plane at ^ ~ 230° and 
then arcs up to i! ~ 140° and b ^ +30°. The best region 
to compute its vertical height is between 140° < £ < 180° 
because it is further away from the Galactic disk than 
elsewhere. We find there that the Northern Arc has a 
Gaussian profile (see Figure 7), centred on 19.9° ± 1.0° 
and with a 2.2 ± 0.4 kpc vertical width. However, high 
extinction regions could be responsible for the decrease of 
star counts under b ~ 15° and the vertical width could be 
underestimated. The thickness of the structure was also 
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Figure 4. As Figure 3, but showing the difference in M-giant 
starcounts between the two hemispheres. The Canis Major over- 
density at {t ~ 240°, b ~ —10°) is clearly visible on the upper two 
panels. The middle panel also shows a fainter population arcing 

over the sky between 110° ^ ^ ^ 210° in the southern hemisphere. 
Another huge arc-like structure is present in the northern hemi- 
sphere in the bottom panel ranging from i ^ 140° to £ ^ 220°. 
In the top, middle and bottom panels, a black pixel corresponds 
to 40, 8 and 8 counts, respectively. 
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Figure 5. As Figure 4, with an overlay of the main asymmetries 
in the Galactic M-giant distribution detected in Figures 3 and 4. 
The strongest of these is the Canis Major over-density, followed 

by the Northern Arc. A fainter Southern Arc is also visible, as are 
two structures at £ = 70° , and £ = 300° , which we label structures 
A and B respectively. 
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Figure 6. Vertical and radial profiles of the Canis Major object, 
derived from sample A. The northern hemisphere is used as a ref- 
erence and its star counts are subtracted from those of the south- 
ern hemisphere. The left-hand panel shows the latitude distribu- 
tion of stars with 15.2 < (m - M)gc < 15.7 and 220° <£ < 260°. 
The lack of data for |6| < 7° puts a high degree of uncertainty 
on the fits, but the best Gaussian fit (full line) is centred on 
b = -6.3° ± 1.0° with a FWHM of 13.1° ± 1.3°, correspond- 
ing to a latitude width of 1.6 ± 0.2 kpc. We also calculated the 
best exponential fit of the population (dashed line): it has a scale 
height of 0.73 ± 0.05 kpc. The right-hand panel shows the radial 
distribution of stars for the region of the centre of the object 
(230° < £ < 250° and -12° < 6 < -7°). The best Gaussian fit 
has been superposed on the star counts: it corresponds to a Gaus- 
sian centred on Dq = 7.1 ± 1.3 kpc, with a 4.2 ±0.3 kpc FWHM. 
We used the Heliocentric distance distribution here because our 
angle selection was made in Heliocentric coordinates, more suited 
for observations along these lines of sight. 



computed, this time in the region where the over-density 
is most pronounced in the northern hemisphere, that is for 
140° < £ < 180° and -hlO° < 6 < +25°. The right-hand 
panel of Figure 7 shows the best Gaussian fit for the radial 
profile of the structure: centred on Doc = 18.1 ± 0.1 kpc, it 
has a FWHM line of sight thickness of 2.2 ± 0.3 kpc. This 
structure is almost certainly due to the M-giants present in 
the "Galactic Ring" of 103 and Y03. Even with the rather 
noisy values derived here, the Gaussian vertical FWHM 
falls in the range of the Y03 and 103 exponential scale 
height measurements (respectively < 4 kpc and ~ 2 kpc) . 

A second southern over-density, the "Southern Arc" in 
Figure 5, ranges over 110° < £ < 210° and 12.0 kpc < 
Dgc < 15.1 kpc (i.e. 15.4 < (m - M)gc < 15.9). This struc- 
ture is detected at S/N = 8. It is possible that its vertical 
distribution has a similar Gaussian profile to that of the 
Northern Arc, but the large highly dust-contaminated re- 
gion near the disk hampers the interpretation of the obser- 
vations. This feature may correspond to the southern hemi- 
sphere population reported by Y03 in the longitude range 
180° <£< 227° and at a distance of Roc = 20± 2 kpc. The 
30% systematic distance error mentioned earlier may be re- 
sponsible for the difference in these two distance estimates. 

The two main over-densities can be followed around half 
of the Milky Way. To diminish the projection problems that 
appear when observing in the quadrant towards the Galac- 
tic centre (i.e. for i > 270° and £ < 90°), we switched to 
the Galactoccntric cylindrical coordinates {R, $, z). We also 
limited our sample to < 30° since the sources are fur- 
ther away from the Sun. The left panel of Figure 8 shows 
the Galactocentric distance modulus distribution of stars 
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Figure 7. Profiles of the Northern Arc structure. The southern 
hemisphere is used as a reference and its star counts are sub- 
tracted from those of the northern hemisphere. The left-hand 
panel shows the latitude profile of the structure for 140° < £ < 
180° and 16.1 < (m - M)gc < 16.7, and the best Gaussian fit of 
the star counts. The distribution has a mean value of 19.9° ±1.0° 
and a FWHM of 12.7° ± 2.3°. At the Hefiocentric distance of 
the object (~ 10 kpc), it corresponds to a vertical width of 
2.2 ±0.4 kpc. The right-hand panel shows the radial profile of the 
structure for 140° <£< 180° and +10° < 6 < +25°. The distri- 
bution is centred on Dec = 18.1 ± O.lkpc with a 2.2 ± O.Skpc 
thickness. Since here the observations are in the direction of the 
anticentre, we were able to use Galactocentric distances. 



for 230° < * < 270° where the prolongation of the Car 
nis Major over-density is visible despite the higher num- 
ber of disk stars. It ranges up to Dec ~ 14.5 kpc from 
at least Dec = 5 kpc but the presence of the bulge pre- 
vents the determination of a precise lower limit. The high 
radial extent of the over-density is certainly due to the angle 
of observation. Indeed, it is possible that the structure is, 
at least partially, not observed perpendicularly but rather 
tangentially. This would also explain its broad profile and 
seems to be confirmed when taking smaller ranges of 
The right-hand panel corresponds to 90° < $ < 130° and 
shows the prolongation of the Northern Arc. It ranges over 
7.2 kpc < Dgc ^ 16.6 kpc which is nearer from the Galactic 
centre than it was in the direction of the anticentre. It is 
surprising to find that it is also denser. This could also be 
due to the angle of observation: if here too we are partially 
observing the structure tangentially, the star counts will be 
summed up along the line of sight. 

To summarize, we were able to detect two substantial 
over-densities of M-giants, the Canis Major over-density and 
the Northern Arc. which together surround half of the Milky 
Way, at varying Galactocentric radius. This population can 
be as close as 5 kpc from the Galactic centre for $ = 90° and 
20 kpc away in the direction of the anticentre. We also find 
it probable that it passes from the northern to the southern 
hemisphere at around £ = 230° . 



3.2 Colour-Magnitude Diagrams 

We now examine the stellar populations of the stars in the 
structure, using the near infrared 2MASS photometry. As 
can be discerned in Figure 9, the (J — Kg, Kg) CMD of a 
5° X 2° region at the centre of the Canis Major object 
shows two striking features: a red-giant branch population 
which does not appear on the reference CMD on the other 
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Figure 8. The Galactocentric distance modulus (m — M)gc dis- 
tribution for M-giant stars with |b| < 30°, E(B - V) < 0.555 and 
for different ranges of the Galactocentric cylindrical coordinate <&. 
The star distribution in the southern hemisphere (full line) is com- 
pared to the northern hemisphere distribution (dashed line); the 
over-densities are also plotted (dotted line). The 230° < <I> < 270° 
region (left panel), which corresponds to f > 270°, shows the 
prolongation of the Canis Major over-density. With a Galacto- 
centric distance modulus of (m — M)gc < 15.8, it is closer to 
the Galactic centre than the over-density we presented on Fig- 
ure 2; the (m — M)gc = 13-5 lower limit is due to the presence of 
the bulge. The over-density has a maximum star count of ~ 300 
(per interval of 0.1 in distance modulus) at (m — M)gc = 15-1 
and slowly decreases with decreasing distance from the Galactic 
centre. The 90° < $ < 130° region (right panel), which corre- 
sponds to £ < 90°, shows the feature labeled as structure "A" 
in Figure 5, which may be the prolongation of the Northern Arc. 
It is also nearer, ranging over 14.3 < (m — M)gc < 16.1 centred 
around (m — M)gc = 15.1, with a profile that is approximately 
symmetric. 



side of the Galactic plane and a high number of stars for 
0.55 < J - Ks < 0.65 and 12 < < 13.5. For comparison, 
we superimpose on Figure 9 the line of the fiducial RGB se- 
quence of the Sgr dwarf, adjusted to 7.1 kpc, the distance of 
the Canis Major structure determined from the analysis of 
Figure 6. Inspection of the subtracted Hess-diagram shown 
in the right-hand panel of Figure 9 reveals that the popula- 
tion in this central field 5° x 2° is slightly (~ 1.5 kpc) more 
distant, than estimated above. An alternative explanation 
could be that there is a stellar population gradient towards 
the centre of the structure. 

This excess RGB population does not appear to be spa- 
tially separated from the Galactic disk, that is, we do not 
detect a magnitude gap separating this population and the 
disk red giants in the data presented in the left-hand panel of 
Figure 9. The structure therefore lies within, or at the edge, 
of the Galactic disk. The other feature that is present at the 
bottom of this red giant branch (at 0.55 < J — Ka < 0.65, 
Kg ~ 12.5) is the corresponding red clump of the stellar 
population. An additional blue plume at J — Kg < 0.5 is 
seen in the subtracted Hess diagram; this feature could be 
an intermediate-age main-sequence population of the Canis 
Major galaxy, although it could possibly also be an artifact 




J-Kj J-Kg J-K_ 

Figure 9. Near infrared colour-magnitude diagrams (from the 
proper-motion cleaned sample B) of the centre of the Canis Major 
object (239° <t< 244° and -10° < fe < -8°) and its symmetric 
region in the northern hemisphere (239° < t < 244° and 8° < 
b < 10°). The Canis Major object CMD shows two features which 
are not in the reference CMD, and that are clearly seen in the 
subtracted Hess diagram shown on the right-hand panel (a black 
pixel corresponds to a count of 25 stars). A red-giant branch is 
present from (J - Ka ~ 0.7, Kg ~ 12) to (J - Ks ~ 1.2, Kg ~ 8), 
while the dense feature at (0.55 < J — Kg < 0.65, Kg ~ 12.5) is 
the red clump. As in Figure 1, the dashed lines show part of the 
colour-magnitude selection region. The solid lines show the RGB 
fiducial of the Sgr dwarf galaxy (left-hand panel of Figure 1), 
adjusted to a distance of 7.1 kpc (the distance to the structure 
determined from the analysis presented in Figure 6). 

of the subtraction. A deeper optical wide-field survey of the 
region may help resolve this issue. 

The colour-magnitude diagram of the Northern Arc re- 
veals a sparse red giant branch, consistent with a Galacto- 
centric distance of ~ 18.1 kpc, as we show in the left-hand 
panel of Figure 10. It appears that there is a small gap be- 
tween the RGB of this population, and the disk RGB stars, 
indicating that the structure is separated from the disk at 
its most distant point. The middle panel of Figure 10 shows 
a comparison region below the Galactic plane; here, we may 
also have detected a similar RGB population, though the 
signal is very weak. This feature may correspond to the 
"Southern Arc" structure highlighted in Figure 5. 

We also searched for the RGB population of the struc- 
ture in small low extinction windows at very low latitude. 
However, the low number of the M-giants of interest com- 
pared to the high number of stars in the disk prevents any 
further conclusion on the connection between the Northern 
Arc and the Canis Major over-densities. The RGB popula- 
tions for 230° < <!> < 270° and 90° < $ < 130° also have 
low contrast because the structures are mixed with the disk 
and the higher number of disk stars weakens the signal due 
to the over-densities. 

3.3 Extinction 

Coul d these inferred over-densities simply be due to errors 
in the lSchlegel. Finkbeiner fc Davij (^9^ extinction maps? 
The colour of the vertical CMD feature of disk dwarfs in the 



A dwarf galaxy remnant in Canis Major 7 




0.5 1 0.5 I 0.5 1 

J - K, J - K, J - K, 



Figure 10. The left-hand panel shows the near infrared colour- 
magnitude diagram (from the proper-motion cleaned sample B) 
of a 15° X 15° region in the Northern Arc over-density (140° < 
t < 155° and 10° < 6 < 25°). Since the Northern Arc is fainter 
than the Canis Major structure, it is necessary to gather sources 
over a wider area to obtain sufficient signal to render the excess 
M-giant population visible. For clarity, we use a Hess diagram, 
except in the colour-magnitude region J — Ks > 0.85, Kg < 12.0. 
The solid line shows the position of the fiducial Sgr dwarf RGB, 
shifted to an equivalent Galactocentric distance of 18.1 kpc, as 
indicated by the analysis presented in Figure 7. A population of 
stars is clumped around this line, indicating the presence of an 
RGB at this distance. A gap between the RGB of the Northern 
Arc over-density and the RGBs of the disk can also be perceived, 
indicating that the structure is separated from the disk in this 
direction. The middle panel shows the Hess diagram and CMD 
of the symmetric region about the Galactic Plane (140° < £ < 
155° and -25° < fe < -10°). A hint of the same RGB is also 
present, though the population is much less numerous; it probably 
corresponds to the "Southern Arc" identified in Figure 5. Due to 
the low contrast of the population over the numerous foreground, 
it proved hard to construct a clear subtracted Hess diagram. Our 
best attempt, shown on the right-hand panel, reveals the presence 
of an RGB-like population beneath the solid line; here a black 
pixel corresponds to 5 stars. The red clump population, however, 
is lost in the noise of the local disk red dwarf sequence. 

extinction-corrected samples is almost identical (to within 
~ 0.1 magnitudes) North and South of the Galactic plane. 
This is demonstrated in Figure 9, where the sequence of disk 
dwarfs is seen both on the left-hand panel at J — Ks ~ 0.6 
in the over-density at {£ = 242°, 6 — —9°), and at the 
same colour in the symmetric field about the Galactic Plane 
{£ = 242°, b = 9°). Additionally, this sequence is detected 
in low extinction regions at t he same J — Ks colour, giving 
us fur ther confidence that thelSchlegel. Fi nkbeiner fc Davis! 
il998h extinction maps are not significantly in error along 
these sight-lines. 

For the RGB-like population present on the subtracted 
Hess diagram on the right-hand panel in Figure 9 to be 
caused by extinction errors, would require a contaminating 
population that runs parallel to that RGB population (but 
both bluer and brighter). The only available population that 
has this property is the RGB of the Galactic disk. However, 
inspection of Figures 2 and 8 shows that the distribution of 
M-giants on opposite sides of the Galactic Plane match well 



at the lower distance limits and at the upper distance limits 
of those diagrams. It is at intermediate distances that a dis- 
crepancy is detected. Extinction errors can only rearrange 
the distance distributions, they cannot cause an apparent 
overabundance of stars all along the RGB. For these reasons 
it is unlikely that the detected over-densities are extinction 
artifacts. 

3.4 Mass of the Canis Major structure 

To estimate the mass of the Canis Major object, we com- 
pare the number of detected M-giants to those observed 
in the Sagittarius dwarf galaxy. For Canis Major, we sub- 
tract the northern hemisphere sample A data from the 
southern hemisphere and sum under the profiles shown in 
Figure 6. Within a 10° radius region (but which misses 
~ 40% of the area due to the |6| < 5° cutoff and due 
to the E(B - V) < 0.555 reddening limit) we detect 2300 
M-giants. The same measurement around the core of the 
Sagittarius dwarf galaxy gives an almost identical number 
of 2200 M-giants. With the caveat that the stellar popula- 
tions are not the same, this suggests that the Canis Ma- 
jor object has a similar total luminosity to the core of the 
Sagittarius dwarf galaxy, which has been measured to be 
Mv = —13.27 (Maiewski et al. 200.3;). If we further assume 
that the mass to light ratio of the two galaxies is simi- 
lar, we conclude that the Canis Major object has an al- 
most identical mass to that of the core of the Sagittarius 
dwarf. The mass o f the latter is a matter of some debate 
ilbata et al. I 119971 llbata fc Lewis | [19 98: H e liTij fc Whit3 
I2OOII: iGomez-Flechoso. Fux fc Martinetlll999^ . but a range 
between 10® M© to 10^ M© seems plausible, with the upper 
bound being the likely pre-disruption mass of that galaxy. 



4 A CLUSTERING OF GLOBULAR CLUSTERS 

While searching for structures in the Galactic globular 
cluster system that could be c onnected with the stream 
of the Sagittarius galax y (iBellazzini. Ferraro fc Ib ata 200^: 
iBellazzini et al.l l2003bll . we noted an interesting group of 
clusters that are not related wit h Sgr and that are re - 
markably confined in phase space (IBellazzini et al.ll2003d) . 
These four clusters, NGC 1851, NGC 1904, NGC 2298 
and NGC 2808, are quite close to one another in space, 
indeed, they may be enclosed in a sphere with radius 
7? ~ 6 kpc centred at Galactic Cartesian coordinates 
{X,Y,Z) = (-11.5, -8.9, -4.5) kpc. Among the outer 
halo globular clusters (those w ith Roc > 10 kpc, see 
IBellazzini. Ferraro fc Ibatal l2003l) there is only one other 
group of globulars that shows similar clustering properties, 
the clusters associated with the main body of the Sgr dSph 
galaxy (M 54, Ter 7, Ter 8, and Arp 2). In Figure 11 the spa- 
tial distributions of the two groups are compared in two dif- 
ferent projections of Cartesian Galactocentric coordinates: 
NGC 1851, NGC 1904, NGC 2298 and NGC 2808 are as 
close one to another as the Sgr globulars. It is interesting 
to recall that the similar proximity of the four Sgr clus- 
ter was not noted until the discovery of their parent galaxy 
jBellazzini et alJl20033) . 

If we consider the radial velocities, the present Sgr 
globular cluster system appears much more compact 
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Figure 11. Comparison of the clustering properties in space 
of the group of globulars NGC 1851, NGC 1904, NGC 2298, 
NGC 2808, NGC 1261 (left panels) and of the globular clusters as- 
sociated with the main body of the Sgr galaxy (right panels; M 54, 
Ter 7, Ter 8, and Arp 2). The clusters are plotted in Cartesian 
Galactocentric projections (upper panels: YX; lower panels: YZ). 
In this coordinate system the Galactic Centre is at (X, Y, Z) = 
(0, 0, 0) kpc and the Sun is at {X, Y, Z) = (-8, 0, 0) kpc. The small 
cross indicates the Galactic Centre and the dotted circles have 
Galactocentric radius from _R = 5 kpc to i? = 30 kpc, in steps of 
5 kpc. Both groups are enclosed in a (continuous) circle having 
radius r = 6 kpc. 

{<tvlsr ~ 21.5 kms~^) than the considered group {ctvi^sr ~ 
95.3 kms~^). However the data presented in Figure 12 show 
that the LSR-corrected radial velocity of the group mem- 
bers correlates well with Galactic longitude. Moreover, their 
Galactocentric distance also correlates with Galactic longi- 
tude. These facts are strongly suggestive of a coherent mo- 
tion along a common orbit. 

After the discovery of this remarkable phase space 
grouping of globular clusters we noted that their positions 
are virtually coincident with the Canis Major structure dis- 
cussed in this paper. The strong analogy with the Sgr case 
brings us to the conclusion that the two structures are as- 
sociated and that NGC 1851, NGC 1904, NGC 2298 and 
NGC 2808 are likely the local remnant of the globular clus- 
ter system of the galaxy that was the progenitor of the Canis 
Major object and the Northern Arc. 

These clusters also show some interesting similar- 
ities in their intrinsic properties. For instance, their 
metallicity is enclosed in the relatively limited range 
-1.03 < [Fe/H] < -1.71, aU of them show an extended 
blue Horizontal Branch in their colour-magnitude dia- 
grams, and NGC 1851 and NGC 2808 are the only two 
halo clusters showing simultaneously a well populat ed red 
clump and a blue tail in their HB morpholo gy (see iHarrisI 
ll996l:lBedin et alJ2000l:lBellazzini et alJ200ll . and references 
therein). 



Figure 12. The upper panel shows the correlation between 
Galactic longitude and the radial velocity (corrected for the mo- 
tion of the Local Standard of Rest) for the clusters associated with 
the Canis Major object. The lower panel presents the correlation 
between Galactic longitude and Galactocentric distance. 

There are other globular clusters that may be related 
to the Canis Major structure. NGC 1261 is also nearby in 
space to the above described group but its velocity does not 
fit well with the velocity gradients shown in Figure 12. We 
will return to this issue in Section 5 below. 

It is also interesting to study the distribution of open 
clusters in light of the discovery of the Canis Major object. 
Figure 13 shows t he open clusters in the WEBDA database 
^Mermi^iodll^95^ ^ with Dq > 6 kpc; a puzzling asymmetry 
is evident with absolutely no clusters for the i < 180° part of 
the Milky Way. Caveats may include strong uncertainties in 
distance and/or age since the directions are often obscured 
and the cluster population is small. Moreover, the presence 
of the Outer Arm of the Galaxy could provide a stronger 
background in the 0° < ^ < 180° quadrant against which 
it could be more difficult to identify clusters. On the other 
hand, part of the 180° < £ < 270° quadrant is behind the 
Perseus Arm which is much closer to the Sun: Therefore the 
considered region suffers from a less severe contamination 
and the identification of remote clusters is easier. 

With these caveats in mind, it is still possible that 
part of these remote open clusters belong to the Canis Ma- 
jor galaxy. In particular Dol 25 {£ = 211.9°, fe = -1.29°) 
and H 18 (£ = 243.1°, fe = 0.44°) which are both young 
(log(age) < 8.5) and with a Heliocentric distance that corre- 
sponds to the Canis Major object could be the consequence 
of a strong perturbation of the disk by the accreted dwarf 
Galaxy and the resulting enhanced star-formation activity 

^ The WEBDA database, developed and main- 
tained by J.-C. Mermilliod, can be found on 
|http: / / obswww.unige.ch/webda/webda.htinT| 
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Figure 13. Distribution of open clusters with Dq > 6 kpc. 
The Galactic centre has been marked (cross) as well as the Sun 
(dotted circle). The full line circles correspond to a distance of 
15 kpc and 20 kpc from the Galactic centre; the dotted circle 
corresponds to the Dq = 6 kpc limit. The symbols correspond 
to the estimated age of the open clusters: log(age) < 8.5 (stars), 
8.5 < log(age) < 9.0 (empty triangles) and log(age) > 9.0 (filled 
triangles). All the open clusters have j6| < 18°. 

in the disk. We also note that the position of the old open 
cluster AM2 (previously classified as a globular) nearly co- 
incides with that of the central peak of the Canis Major 
structure, within the (large) distance uncertainties. Unfor- 
tunately the lack of any radial velocity estimate for this 
cluster prevents any definite conclusion about its possible 
physical association with the newly discovered object. 



5 N-BODY SIMULATION 

We undertook a series of N-body experiments in an attempt 
to determine whether a single accretion of a dwarf galaxy 
could be responsible for the detected over-densities and the 
correlations seen among globular and open clusters. In this 
first study, we concentrate on constructing a simple model 
to guide interpretation of the observations, and will defer a 
more detailed analysis to a subsequent contribution. 

For our simulations, we take a static Milky Way poten- 
tial, which is not perturbed by the passage of the accreted 
dwarf g alaxy. The potential model 2b of|Pehnen & Binnev 
Jl998ai) is adopted, but with the halo component modeled 
by a spherical NFW halo. Taking a static potential is unre- 
alistic, of course, but it saves a great deal of computation 
time, making the simulation tractable. The key effects that 
this approach leaves out are the back-reaction of the dwarf 
galaxy on the Galactic disk, and the back-reaction on the 
halo. The disk will likely become deformed (though note that 
the progenitor of the Galactic "Ring" has an angular mo- 
mentum vector that is almost at right angles to that of the 



Galactic warp, so it is unlikely to warp the disk strongly), 
and these deformations will in turn alter the trajectory of 
the accreted object. Since we do not include a live halo com- 
ponent, we neglect the back-reaction on the halo (formation 
of a wake), and the resulting dynamical friction on the dwarf 
galaxy. 

A useful working picture is that the dwarf galaxy began 
its life as a typical dark matter-dominated dwarf satellite 
galaxy, far away in the Milky Way halo. The dark mat- 
ter halo of the dwarf galaxy became progressively more and 
more disrupted as the orbit decayed. However, the central 
stellar component stayed largely intact during this initial 
phase. We start our simulation 2 Gyr ago, and assume that 
by this stage much of the dark matter has been removed 
from the system, leaving a much less massive dwarf galaxy 
that will not suffer significant dynamical friction. We there- 
fore do not add dynamical friction in our simulation. 

The simulations were performed with falcQN, a fas t 
and momentum-conserving tree code (lDehnenll2nnrt |22q3) ■ 
Gravity was softened with the kernel "P2", with a soften- 
ing length of 0.1 kpc. The minimum time-step was set to 
2-14 _ g ^ 10~^ Gyr, and the tolerance parameter 6 — 0.6. 
In order to obtain a starting position and velocity for the 
dwarf galaxy model, we refined iteratively an initial guess 
of these values, and required that the resulting orbit fit the 
velocity and position information of the Y03 SDSS and the 
2MASS fields presented above, as well as the globular clus- 
ters. Acceptable orbits could be fit for both a co-rotating 
and a counter-rotating object, so we examine both cases be- 
low. These fitted orbits have an azimuthal period of 0.45 Gyr 
and 0.3 Gyr, respectively, for the co-rotating and counter- 
rotating models. 

The initial dwarf galaxy was modeled with a King model 
with 10^ particles. We experimented with several such mod- 
els, trying out mas ses in the range o f 2 x 10* to 10^ Mq (the 
range suggested bv llbata et alJ200,?l . and which is consistent 
with the mass estimate derived above), and with different 
tidal radii and concentration parameters. Some fine-tuning 
of these parameters is needed to give both a stream and 
remnant that approximates the observed distribution. One 
of the better fitting models had a mass 5 x 10* M0, tidal 
radius 2.5 kpc and WQ = 4.5, and it is this model that we 
will discuss below. However, these values are by no means 
the required initial parameters of the satellite galaxy. They 
are sensitive to the orbit, and must therefore also be sen- 
sitive to the decay rate of the satellite due to dynamical 
friction, which we have neglected. The model is therefore 
only illustrative. 

In the top panels of Figures 14 and 15, we show respec- 
tively the x-y configuration (looking down onto the Milky 
Way from the north Galactic pole) at the end of the simula- 
tion, for the co-rotating and counter-rotating models. Given 
the simplicity of the assumptions the similarity between the 
simulations and observations is remarkable, though we feel 
that a formal statistical comparison is not very meaning- 
ful at this stage due to the fact that we have left out the 
self-gravity of the Milky Way. However, what the simula- 
tions do show is that the streams that arise naturally dur- 
ing the tidal disruption of a dwarf galaxy can be pulled 
out to be seen as multiple streams, even in the relatively 
short time (~ 2 Gyr) that we simulated. This may account 
for the different peaks in the distance modulus distribu- 
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tions (Figures 2 and 8). Figures 14 and 15 also show the 
positions of the globular clusters that have Galactocentric 
distances Roc > 8 kpc, and that lie within 8 kpc of the 
Galactic plane. To select the closest star clusters to our sim- 
ulation, we choose somewhat arbitrarily a simple statistic: 

^2 ^ |g-g^odeih ^ 



; we list the globular 



2 kpc ; ' V 20kms-l 

and open clusters that have x < 3 in the prograde and ret- 
rograde models, in Tables 1 and 2 respectively. We defer 
analyzing the likelihood of the alignment of these star clus- 
ters with the Canis Major stream to a future contribution 
when we present a more realistic simulation with a live Milky 
Way. 

The distribution of the simulation particles projected 
onto the sky is shown in Figure 16; the models nicely repro- 
duce the curving behaviour of the Northern Arc observed 
towards the Galactic anticentre (compare to Figure 4), the 
Canis Major over-density as well as the over-densities "A" 
and "B" (see Figure 5) seen towards £ = 70° and e = 300°. 
The Heliocentric radial velocity profiles of the models are 
given in Figure 17. To our surprize, we find that current 
data do not appear to prefer the prograde model over the 
retrograde. Future measurements of the radial velocity of the 
stream in directions sufficiently far away from the Galactic 
anticentre should resolve this uncertainty. 

The numerical simulations presented here lend support 
to the notion that the Canis Major over-density is the rem- 
nant of a larger accreted galaxy whose tidal disruption gave 
rise to the asymmetries in the Galactic M-giant distribu- 
tions that we report above. If we accept this model, we find 
that the progenitor had an orbit that was closely aligned 
with the Galactic Plane (the orbit is inclined at only ~ 10° 
from the Plane). It is unclear at present whether the object 
follows a prograde or a retrograde orbit, though a-priori the 
prograde case would appear to be more likely since it would 
allow for stronger interactions with the Galactic disk, due 
to which the orbit would have decayed to its present, rather 
circular, state. The remaining debris has an apocentre of 
^ 20 kpc and a pericentre close to the Solar radius. Indeed, 
according to the simulations displayed in Figures 14 and 15, 
the debris from the Canis Major dwarf streams through the 
Solar Neighbourhood. Comparing the kinemati cs of the Ca- 
nis M ajor stream to the stream discovered bv lHelmi et alJ 
1^9^, the only other large stream known to pass close to 
the Sun, shows that the streams of our models have very 
different angular momentum. The angular momentum dis- 
tribution of the particles within 2 kpc of the Solar Neigh- 
bourhood is shown in Figure 18; these particles display very 
different kin e matic s, and are unlikely to be associated to the 
iHelmi et al.l il999ll stream. 



6 DISCUSSION AND CONCLUSIONS 

We have shown that the distribution of M-giant stars at 
low latitude is highly asymmetric. The most striking over- 
density is due to an elliptical object present in the disk, in 
the direction of the Canis Major constellation in the south- 
ern Galactic hemisphere at Dec = 12.0 ± 1.2 kpc. This 
structure contains as many M-giant stars as the remnant 
of the Sagittarius dwarf galaxy. A less dense structure, but 
of huge proportions, the Northern Arc curves around the 
Galaxy in the direction of the Galactic anticentre. This arc 



"1 1 1 1 t;^ 

Prograde . 




1 1 ^ 

-20 


_^ — 1 1 - . ^rfffmam^ 

-10 


10 


20 




K (kpc) 






^^^^^^^^^^^^^^ 


-ZDD 


-tat a 


1DD 


sua 







Figure 14. Two different projections of an N-body simulation 
of the disruption of the progenitor of the Galactic "Ring" and 
the Canis Major structure. The Galactic Centre is at the ori- 
gin, and the Sun at {x,y,z = —8,0,0 kpc). The dwarf galaxy 
model was placed on an orbit that co-rotates with the Galactic 
disk. The upper panel shows the x-y distribution (i.e. as seen 
from the Galactic north pole), with colour-coded radial veloc- 
ities. In black, we show the positions and distances of the INT 
fields from Ibata et al. 1 2003) (large open circles), and the 2MASS 
data presented in this contribution (open squares and diamonds). 
The diamonds represent those fields for which a distance esti- 
mate was secured by summing in Galactocentric, rather than 
Heliocentric distance. The colour-coded graph markers show the 
positions, di s tances and velocities of the four SDSS fields from 
lYannv et alj j20o3) (open triangles), and of the 4 globular clus- 
ters detailed in section 4 (star symbols). The filled circles rep- 
resent all globular and open clusters of known radial velocity 
within 8 kpc of the Galactic plane, and that have Galactocen- 
tric distances Rqc > 8 kpc; the subset of these that are circled 
black are those objects listed in Table 1, which are close in phase- 
space to the simulation particles. This velocity-distance-position 
coincidence with the model suggests that several globular clusters 
around the Milky Way may be associated with the Canis Major 
accretion event. 



A dwarf galaxy remnant in Canis Major 1 1 




1 1 

-20 


1 1 1 1 1 

-10 10 


1 1 

20 




K (kpc) 






Hi ' 1 ' 1 

^^Kj .1.1 





2DD WD 



Figure 15. As Figure 14, but for a model on a retrograde orbit. 





Figure 17. The Heliocentric radial velocity of the prograde 
model (upper panel) and retrograde model (lower panel) as a 
function of Galactic longitude, is compared to the Y03 SDSS fields 
(open triangles) and the 4 globular clusters NGC 1851, NGC 1904, 
NGC 2298, and NGC 2808 (star symbols). In addition, the open 
circles show, from right to left, the positions of the star clusters 
in Table 1 (for the top panel) and Table 2 (for the bottom panel) , 
in the order listed. 




2000 
J (kpc km/s) 



4000 



Figure 16. The £—6 distribution of the prograde model (upper 
panel) and retrograde model (lower panel). Most of the observed 
large-scale asymmetries of the Galactic M-giant distribution are 
reproduced in this projection of the simulations: the broad curved 
distribution of the Northern Arc seen towards the anticentre, the 
strong over-density seen towards Canis Major at £ ~ 240°, 6 ~ 
— 10° , a weak feature similar to the "Southern Arc" is also present, 
plus the over-densities "A" and "B" at <I> ~ 100° , and <E> ~ 260° . 



Figure 18. The larger dots show the angular momentum dis- 
tribution of particles in the prograde simulation within 2 kpc of 
the Sun, while the smaller dots show the same parameters for 
particles in the retrograde simul ation. These distributions differ 
substantially from that found bv lHelmi et alJ jl999t) for the star 
stream detected in the Solar Neighbourhood (star symbols). The 
particles in both the Canis Major simulations follow a more pla- 
nar orbit, with the prograde particles having typical Galactic disk 
kinematics. This large difference implies that the C anis Major 
stream is not associated with the lHelmi et al.l il999l) structure. 
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ranges from the Galactic plane ai £ ^ 230° to b ~ +30° at 
e ~ 140°. Between 140° < £ < 180° it has a mean displace- 
ment of ~ 20° above the Galactic Plane, with a distance of 
Dec ~ 18.1 ± 0.1 kpc, it is separated from the disk in the 
radial direction along the line of sight, and has a Gaussian 
profile of width 11.6° ± 1.9°. 

Several other structures are also detected. In the south- 
ern hemisphere, the "Southern Arc" population with Galac- 
tocentric distance 12.0 kpc < Dqc < 15.1 kpc is observed 
between 110° < £ < 210°; this may be the continuation 
of the southern hemisphere structure reported by Y03 for 
Rgc = 20 ± 2 kpc and 180° < £ < 227°. Two other over- 
densities are detected towards the Galactic bulge, and close 
to the Galactic Plane: one is found between 230° < $ < 270° 
at Dec < 14.5 kpc, while the other is between 90° < $ < 
130° at Dae = 10.5 kpc. 

Given these observations, we examine the idea of an 
in-plane accretion of a dwarf satellite galaxy onto the 
Milky Way. The realization that the four globular clusters 
NGC 1851, NGC 1904, NGC 2298, and NGC 2808, are close 
to each other in phase space and are also close in posi- 
tion and distance to the Canis Major over-density also lends 
strong support to this possibility. In this scenario the rem- 
nant of the satellite is likely the object in the direction of 
Canis Major. During the tidal disruption process, the dwarf 
galaxy loses stars in a long stream which loops around the 
Milky Way and is now seen as the Monoceros "Ring" struc- 
ture. The two other hypotheses presented by 103 do not seem 
to be allowed by our observations. Indeed, the height of the 
structure (generally more than 3 kpc) could hardly be ex- 
plained by a spiral arm. Furthermore, the multiple peaks 
in the distance distribution are more suggestive of multi- 
ple crossings of a stream rather than of a ring-like Galactic 
structure, produced by the dynamical evolution of ancient 
stellar warps. 

The simple N-body models we presented reproduce 
many of the observed spatial features, including the Ca- 
nis Major over-density, the vast ring-like arc towards the 
Galactic anticentre, the over-densities towards £ — 70° and 
£ = 290°, as well as faint debris scattered out to large dis- 
tance below the Galactic Plane in the anticentre direction. 
The models also reproduce the available kinematic data and 
point to the possibility that several star clusters, apart from 
those close to the Canis Major over-density, could be re- 
lated to the Canis Major Stream. The success of this simple 
model strongly suggests that many, if not all, of the asym- 
metries in the Galactic M-giant population presented above, 
are due to a single accretion event. The accreted satellite 
must have had an orbital plane closely aligned to the plane 
of the Mi lky Way galaxy. Ou r analysis shows that the "shell" 
model of lHelmi et all J2003l) gives poor agreement with the 
observations, but we confirm the tidal stream model for the 
"Monoceros Ring" proposed by Y03. 

The observations and model also suggest a connection 
to the Galactic thick disk. The orbital eccentricity of the 
model is e = {rapo - rperi)/{rapo + Vperi) ~ 0.5, a value 
that is typical for Galactic thick disk stars of metallicity 
[Fe/H] ~ -1.2 Jchib a & Beera,.200a) . Furthermore, the ver- 
tical exponential scale height of the debris 0.73 ± 0.05 kpc, 
is very close to t he th ick disk value of 0.8 kpc determined by 
iRevle fc RobinI (1200111 . It is therefore possible that we have 
uncovered one of the galactic building blocks from which the 



thick disk grows. This analysis suggests that the thick disk 
is continually growing, even up to the present time, through 
accretions of dwarf galaxies in co-planar orbits. Although 
with a progenitor mass of ~ 10® M©, as suggested by our 
analysis, only a handful of accretions of the magnitude of 
that of the Canis Major dwarf galaxy are needed to popu- 
late the thick disk entirely. Globular clusters with disk-like 
kinematics may also have been brought to their present lo- 
cations in this way. The observations a nd modeling we pre - 
sented strongly support the analysis of lAbadi et al.l i2002l) . 
in which a large proportion of thick stars have their origin in 
accreted galaxy fragments. If the scenario we have developed 
here is correct, thick disks should be a generic component of 
giant spiral galaxies, since accretions are a generic feature 
of galaxy formation. Furthermore, at the edge of the disk 
where the dynamical times are longest, their structure will 
likely be a complex mess where streaming debris from the 
most recent accretions piles up at apocentre. These messy 
structures may in deed already have been detected in the An- 
dromeda galaxy dlbata et alJl2001bl : lFerguson et al.l Eo02). 

Further observational work is needed to constrain the 
kinematics of the structures we report here. These con- 
straints are needed to determine whether the structure fol- 
lows a prograde or a retrograde orbit, to verify the kine- 
matic predictions of our models, as well as to guide future 
simulations. Detailed abundance distributions will also be 
invaluable in order to study the chemical difference between 
an accreted population, such as that of the Canis Major 
galaxy, and that of the normal Galactic disk. Future simu- 
lations with a live Milky Way will provide a more realistic 
comparison to the observations. It will be interesting to in- 
vestigate the reaction of the Galactic disk stars and in par- 
ticular the reaction of the HI gas, to this sizable accretion, 
while the back-reaction onto the Canis Major dwarf galaxy 
will be essential in order to model its dynamical evolution 
accurately. 
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